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ABSTRACT 

Using numerical simulations at moderate magnetic Reynolds numbers up to 220 it is shown 
that in the kinematic regime, isotropic helical turbulence leads to an alpha effect and a turbu- 
lent diffusivity whose values are independent of the magnetic Reynolds number, Rm, provided 
i?,„ exceeds unity. These turbulent coefficients are also consistent with expectations from the 
first order smoothing approximation. For small values of R-^, alpha and turbulent diffusivity 
are proportional to _/?,„• Over finite time intervals meaningful values of alpha and turbulent 
diffusivity can be obtained even when there is small-scale dynamo action that produces strong 
magnetic fluctuations. This suggests that small-scale dynamo-generated fields do not make a 
correlated contribution to the mean electromotive force. 
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1 INTRODUCTION 

The generation and maintenance of large-scale magnetic fields 
in stars and galaxies is often studied w ithin the framework of 
the m e an-field dynamo (MFD) ; see, e.g., iMoffattI ( Il978h : [Parked 
( ll979h ; lKrause & Radleij ( Il98(]|) . A particularly important driver of 
MFDs is the a-effect. For isotropic turbulence and weak magnetic 
fields, i.e. in the kinematic regime, the a-effect can be expressed 
purely in terms of the kinetic helicity. Research in recent years has 
mostly been concerned with clarifying the effects of nonlinearity, 
but there are serious uncertainties even in the linear (kinematic) 
regime. In particular, whether or not a can then be expressed in 
terms of the kinetic helicity depends on the applicability of the first 
order smoothing approximation (FOSA) or other closures used to 
calculate a. Such approaches become questionable when the mag- 
netic Reynolds number, i?m, is large, i.e. when the magnetic diffu- 
sion time is long compared with the turnover time which, in turn, is 
comparable with the correlation time of the turbulence. In the high 
conductivity limit FOSA can only be applied if the correlation time 
for the velocity field is much smaller than the eddy turnover time. 
This is not the case for high Reynolds number turbulence, where 
the two time scales are equal, i.e. the Str ouhal number is unity 
JBrandenburg & Subramanianl l2005bl |2007|) . so FOSA should in 
principle break down. I n this case all hi gher order terms need to 
be taken into account JKnoblochl 1 197q) . Furthermore, high 7?m 
random flows typically lead to a fluctuation dynamo which leads 
to rapidly growing small-scale magnetic fields independent of the 
mean field. This also breaks the assumption made by FOSA that 
fluctuating fields are much smaller than the mean field. 

The existence of a-effect and turbulent diffusion has been 
worrying dynamo researchers over several decades. Assuming 
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steady flow patterns. Ichildressi ( 1 1 9791) found that motions that con- 
centrate magnetic fields into thin flux sheets lead to an a-effect 

— 1/2 

whose value diminishes with _Rm like _Rm . But for an analo- 
gous helical motion which concentrates the field into an axial flux 
rope, a tends to a finite limit as Rm -^ oo- He conjectured that 
the latter estimate may be typical of steady three-dimensional mo- 
tion. The validity of turbulent dif fusion has been q uestioned by 
IPiddingtonI ( Il98ll) . Calculations bv lKraichnaril ( ll976l) suggest that 
a and turbulent diffusion converge to finite values for statistically 
isotropic veloci ty fields with gaussian statist ics. However, numeri- 
cal simulations I IDrummond & Horganlll986i) using a frozen veloc- 
ity field suggest that in the limit of large magnetic Reynolds num- 
bers a tends to zero. Based on specific imposed (kinematic) flow 
patterns it has been suggested that there is no si mple relation be- 
tween a and helicity of the flow; see lCourvoisier. Hughes & Tobiag 
(2006). In fact, their results may suggest that in the kinematic 
regime, a exhibits a strong Rm dependence (for large Rm up 
to 2 X 10^) and could change sign for _Rm ~ 20 and long 
correlation times. It is important to emphasize that in the non- 
linear regime, i.e. for finite magnetic field strength and includ- 
ing the Lorentz force, a strong R^ dependence is now indee d 
we ll established fcattaneo & Hughe j H 9961; lBrandenbur3 1200 ll) ; 
see lBrandenburg & Subramanian ( 2005a|) for a review. However, in 
the following we shall be concerned with the purely linear regime. 
In order to clarify the i?m dependence in the kinematic regime 
we perform numerical turbulence experiments where we adopt an 
externally imposed body force to drive the flow. This is a com- 
mon technique applied in simulations and helps develop homoge- 
neous and isotropic turbulence which is easier to handle analyti- 
cally. While the calculation of a from a turbulence simulation is 
relatively straightforward by imposing a uniform magnetic field, 
the calculation of turbulent diffusion is more uncertain. One possi- 
bility is to determine the decay rate of an initial l arge-scale mag- 
netic field JYousef. Brandenburg & RUdigenuOOBO . Another more 
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reliable method is to calculate both a and turbulent diffusion ten- 
sors simultaneously by computing the mean electromotive force us- 
ing a number of different nonuniform test fields in different direc - 
tions and of different spatial structure JSchrinner et alJlOOSLuOOTT) . 
This is also the approach used in the present paper. 



2 TEST-FIELD PROCEDURE 

In MFD theory, one averages the induction equation to obtain the 
standard dynamo equation for the mean field B, 

^ = V x(U xB + S-'qV xB), V-B = 0. (1) 

ot 

This averaged equation now has a new term, the mean electromo- 
tive force (emf) £ — u x b, which crucially depends on the sta- 
tistical properties of the small-scale velocity and magnetic fields, 
u and b, respectively. A central closure problem in MFD theories 
is to compute the mean emf £ and express it in terms of the mean 
field itself. Assuming that the mean field is spatially smooth, the 
mean emf £ can then be expressed in terms of the mean magnetic 
field and its first derivative in a manner 



c-i — OiijjDj -\- Tfijj^Iljjf^^ 



(2) 



where atj and rjijk are turbulent transport coefficients written in 
tensorial form, and a comma denotes partial differentiation. 

In numerical simulations, the full aij and rjijk tensors are de- 
termined by first computing the mean emf using t est fields in dif- 
ferent directions and of different spatial structure I Schrinner et al.l 
[2005, 2007). In what follows, we employ xy averages such that 
the resulting mean fields are expressible only as functions of z 
and t. Some details of the tes t field method applie d to this case 
have already been described bv iBrandenburg J2005I) . Note that the 
solenoidality condition then gives, Bz = const — 0. Hence, one 
only needs to compute the four components of atj and rjijs with 
i,j = 1,2. Here the numbers 1, 2, and 3 refer to Cartesian coordi- 
nate directions x, y, z. In order to obtain the 4 + 4 = 8 unknown 
coefficients, we need the x and y components of 4 different test 
fields, Bj , where i, p, and q take values 1 and 2. 

In order to compute the components qh and 77123, for exam- 
ple, it suffices to consider the following two test fields. 



(3) 



Here fci is the smallest non-vanishing wavenumber in the domain, 
and p — 1 or p = 2 denotes whether we take a cosine or sine 
behaviour for the test field, and q = 1 or g = 2 depending on 
whether the non-zero component of Bi is the x-component or y- 
component respectively. We insert these test fields into the relation 








(4) 



Since q — 1, only the j = 1 (x-) component contributes to the sum 
over j above. For the two values of p and with index i unspecified, 
we have 



£i = Qii cos fci2 — 77^13/01 sin fciz, 
£i = an sin kiz -\~ rjii^ki cos k\z. 



(5) 
(6) 



Similarly, for test fields which only have non-zero y-components, 
that is with q = 2, one obtains a similar pair of equations with 
the same arrangement of cosine and sine functions, but with 0^2, 
»7i23 in place of Qa and 77^13- So, for each value of i one obtains 



independent matrix equations for the unknown coefficients aij and 

ilij-j. as 



Jjij'ik-i 
where 

M 



M-f| 



coskiz —smkiz 
sm k\z cosfciz 



(7) 



(8) 



is the same matrix for each value of q and each of the two com- 
ponents i = 1, 2 of £\'^ . Note that det M = 1, so the inversion 
procedure is well behaved and trivial. 

Given the form of the test fields, we can compute the mean 
electromotive force f "^ = u x f"^ for a given test field B . The 
mean emf £'"' is computed by solving the equations 

dbP" 



dt 



= V X C7 X 6*"' + u X B"" + G"" + r,V^b 



2yvq 



(9) 



for each test field B along with the momentum equation for the 
fluctuating velocity field (see further below)|jHere 



G'"' = V xiuxb""" -ux bP") 



(10) 



is a nonlinear term that would normally be neglected under FOSA, 
but will be retained in the numerical simulations. For sufficiently 
large values of /?,„ the small-scale field b'"' can grow exponentially 
due to small-scale dynamo action. An important question to ask is 
whether the predictions of FOSA work even in the presence of such 
a small-scale dynamo. 

We adopt an isothermal equation of state with constant speed 
of sound, Cs, so the momentum and continuity equations are 



du 
'dt 
d\np 



= -u- Vu-c^Vlnp + Z + p ^V-2puS, 



= — u ■ Vlnp — V ■ u, 



(11) 



(12) 



where / is a random forcing function consisting of circularly po- 
larized plane waves with positive helicity and random direction, S 
is the traceless rate of strain tensor The length of the wavevector 
of the forcing function, |fcf |, is chosen to be between 4.5 and 5.5, 
so the average is around fcf — \k{\ ~ 5. The ratio fcf/fci is re- 
ferred to as the scale separation ratio. It must be large enough to 
ensure that higher derivatives in equation ^ can be ignored and 
that a large-scale field could grow, if it was allowed to do so. For 
fully helical turbulen ce, a ratio of 5 is already sufficient, b ut 2.3 
is not; see Fig. 23 of [Haugen. Brandenburg. & Dobleii ( 120041) . We 
adjust the strength of the forcing such that the flow remains clearly 
subsonic (mean Mach number is below 0.2), so for all practical pur- 
poses the flow can be considered nearly incompressible. The details 
of the forcin g function used in the present work can be found in Ap- 
pendix A of JBrandenburg & Subramanianll2005a) . 

We do not include the Lorentz force in the momentum equa- 
tion since we want to study the mean emf in the purely kinematic 
limit. We ignore here the possibility of a mean flow; such flows 
have not been seen to emerge under the simple conditions con- 
sidered here. In the following we use the PENCIL Cooqj, where 
the test field algorithm has already been implemented. We employ 
periodic boundary conditions and use a resolution of up to 512'' 



^ Note that in the corresponding expression (27) of Brandenburg (2005) 
the U term is incorrect. This did not affect his results because U = 0. 

^ |http : //www.nordita. org/ software/pencil- code| 
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Figure 1. Parameters space covered by the simulations presented in this 
paper. 
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Figure 2. Time series of the z averaged values a and r]t for i?m = 220 
and Re = 2.2. Time is expressed in turnover times, t = (wrmsfcf)"^. 
Dash-dotted lines give the time average and dotted lines mark eiTor margins 
obtained by averaging over subsections of the full time series (see text). In 
this run there is small-scale dynamo action leading to strong fluctuations 
for t/r > 25, making the determinations of reliable averages harder at late 
times. 



meshpoints for the run with the largest fluid Reynolds number. We 
use a domain of size {2-k)''^, so the smallest wavenumber is fci = 1. 



3 RESULTS 

We are particularly interested in the dependence of Oij and riij^ 
on Rm and have considered cases where either the fluid Reynolds 
number was fixed, Re = Urms/ii^k{) = 2.2, or where Re = 
10i?ni, corresponding to a magnetic Prandtl number of Pni = 
v/rj = 0.1; see Fig.[T] 

The flow is isotropic and, not surprisingly, we find that, to a 
good approximation, the a and r]t tensors are isotropic with 



ail = Ce22 



?7l23 = — ?7213 = Vt, 



(13) 



and ai2 — 021 = r;ii3 = 77223 — 0. The quantity rjt is simply 
referred to as turbulent magnetic diffusivity. Although the code is 
capable of solving for the full Oij and rjijs tensors, we simplify 
matters by solving only for an and rjns using just two test fields. 
We present the results for a and rjt normalized to the respec- 
tive expressions obtained using FOSA for large magnetic Reynolds 
numbers ( Moffatt 1978 : Krause & Radlei 198Q) , 



QO = 



77to = i^w^, 



(14) 



where r denotes the correlation time of the turbulence. Using the 
definitions of the Strouhal number. 



ot — Tt/rnis"jf 



(15) 



and t he fact that St sa 1 fo r large enough magnetic Reynolds num- 
bers dBrandenburg & Subramanian 2005b . 2007) . we expect 



ao 



rjto = 3M1 



rmsn'f 



(16) 



for a flow that is maximally helical and has positive helicity. For 
Rm < 1 the relevant value of r is no longer the dynamical time 
scale, (urmsfcf)^^, but the resistive one, {rjtk'f)^^. Therefore, both 
a/ao and -qt/rito have to be scaled by i?m for Rm < 1. 

The test field procedure yields a and rit as functions of z and 
t. Since the turbulence is homogeneous, we average these data first 
over 2 and calculate then time averages over the full time series. 



An example of such a time series is shown in Fig.|2]for Re = 2.2 
and Rm ~ 220. Note that even for Rm S> 1 the time averages 
of a/ao and rjt/rito are close to unity, i.e. the predictions from 
FOSA appear to be reasonably accurate. We use such time series 
to calculate error bars as the maximum departure between these 
averages and the averages obtained from one of three equally long 
subsections of the full time series. 

The degree of fluctuations in the time series of a and rjt is 
quite moderate and not at all as strong as in the nonlinear regime 
where fluctuations of a c an even dominate over the mean value 
JCattaneo & Hugheslll996l) . The latter is likely to be a consequence 
of a very small (catastrophically quenched) mean value in the non- 
linear regime. On the other hand, in the kinematic regime the 
level of fluctuations is found not to vary significantly with Rm', 
see Fig. [3] A wea k Rm dependence has a lso been found in the 
presence of shear I IBrandenburg et al.ll2007l) . where such fluctua- 
tions can contribute to dynamo action by an incoherent a effect 
jVishniac & Brandenburg 119971) . However, at late times, i.e. to- 
ward the end of the simulation, the degree of fluctuations increases 
(Fig.|2j. This has to do with the emergence of small-scale dynamo 
action that leads to the production of strong small-scale magnetic 
fields, b'"', although this does not affect the resulting time averaged 
emf for a reasonably long stretch of time. 

In our simulations with Re = 2.2, small-scale dynamo action 
with exponential growth of the rms value of b^'' is found when 
Rm is larger than a certain critical value Rm.a that seems to be 
somewhere between 28 and 44; see Fig. |4] which shows the growth 
of the rms value of b^' for Rm = 220 and Pm = Rm/Re = 100. 

On the other hand, for Rm < Rm,cT the rms value of 
the small-scale field settles to a constant value. In the super- 
critical case the magnetic field is highly intermittent in the 
sense that only in a few places the magnetic field reaches 
large positive and negative values; see Fig. |5] Such inter- 
mittency is typical of large-Pm small-scale dynamo action 
in the kinematic stage iZeldovich. Ruzmaikin & SokoloflJll99(]| ; 



[Brandenburg & Subramanianll2005ah 



The corresponding normalized values of a and rit as a function 
of Rm are shown in Fig. (6) for the case when Re — 2.2. Fig. |7] 
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Figure 3. Dependences of arms (normalized by rjtofei, solid line), and 
Vt,TinB (normalized by rjto, dashed line), on the magnetic Reynolds number 
for Re = 2.2. The vertical bars denote the error estimated by averaging 
over subsections of the full time series. 
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Figure 5. Colour-coded (or grey-scale) representation of 6^^ on the periph- 
ery of the box for a run with Re = 2.2 and R^ = 220 at t/r = 28. The 
colour/grey scale has been clipped at ±100, even though the extrema are 
at ±200. Note the extreme intermittency as evidenced by the presence of 
extended nearly field-free regions. 
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Figure 4. Evolution of the root-mean- squared value of the small-scale mag- 
netic fields b" for Re = 2.2 and fim = 220. (The result for b^^ is over- 
plotted as a dashed line, but it is almost indistinguishable from the solid line 
for b .) Note the nearly exponential growth for t/r > 10, as illustrated 
by the straight line with a slope corresponding to a growth rate of 1.55t~^. 



shows the results for simulations with Pm = 0.1. It turns out that 
in all cases with Rm > 1, ct/ao ~ rjt/rjto ~ 1, while for Rm < 1 
these ratios are equal to Rm. 

Also in the case with Pm ~ 0.1, where we explore larger 
values of Re up to 300, a/ao and rjt/rito reach values close to 
unity. 

We may conclude that for isotropic homogeneous turbulence 
the high conductivity results obtained under FOSA are reasonably 
accurate up to the moderate values of Rm that we have tested. Inter- 
estingly, this conclusion is obtained even in the presence of a small- 
scale dynamo, where b is growing exponentially. Conversely, when 
different values of a and ryt are found under specific circumstances, 
then this must be related to the nature of these circumstances and 
does not indicate the break down of FOSA in general. 

How is it then possible that the high conductivity limit of 
FOSA works even though the correlation time of the velocity field 
is comparable to the eddy turnover time? A possible reason for this 
could be that in the kinematic regime the high conductivity limit 
of FOSA gives similar predictions than the minimal r approxima- 
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Figure 6. Dependence of the normalized values of a and rjt on Rm for 
Re = 2.2. The vertical bars denote twice the error estimated by averaging 
over subsections of the full time series (see text). The run with Rm = 220 
(Re = 2.2) was done at a resolution of 512^ meshpoints. 



tion (MTA) where the triple correlations are not neglected, but re- 
placed by the quadratic correlations divided by a turnover time. 
This closure assumption is not well justified, although numerical 
simulations (for Rm < 300) support some aspects of this closure 
JBrandenburg & SubramanianI l2005bl , llOOTD . Let us also empha- 
size that in this work we have not probed two particular aspects 
where FOSA and MTA depart from each other: feedback from 
helical mean fields (see Brandenburg & Subramanian 2005a, for 
a review), and strong time dependence. In the latter case a time 



© 2006 RAS, MNRAS 000,[T]l5] 



Kinematic alpha effect 



1.00 



0.10 



0.01 




0.01 0.10 



1.00 10.00 



1.00 



0.10 



0.01 




0.01 



Figure 7. Same as Fig.|6] but for P^ = 0.1. The run with Rm ^ 22 
(Re ss 220) required a resolution of 512"^ meshpoints. 



derivative of £ would be important in gi ving the mean field equa- 
tion the character of a wave equation |Blackman & Field 120021 : 
iBrandenburg. Kapvla. & Mohammedll2004) . 



4 CONCLUSIONS 

The present work has shown that for isotropic turbulence the 
first order smoothing results give quite accurate expressions for 
the Q effect and the turbulent magnetic diffusivity in the kine- 
matic regime. This result comes almost as a surprise given that 
in recent years mean field theory has been seriously challenged 
based on numerical simulations. However, we can now clearly 
say that for moderate Reynolds numbers up to about 220 and un- 
der the conditions stated (scale separation, isotropy, etc.) there is 
no evidence that the kinematic results obtained using FOSA are 
flawed, even though its applicability can then no longer be guar- 
anteed and we k now of its shortcomings in the nonlinear regime 
JBrandenburg & Subramanian 2005a ). As explained above, a pos- 
sible reason for this might be that the predictions of FOSA and 
MTA are rather similar, even though FOSA loses its justification 
whilst MTA hinges on a not well justified closure hypothesis. 

The emergence of small-scale dynamo action for supercritical 
values of Rm of about 30 does not seem to affect the average values 
of a and rjt . This suggests that the exponentially growing part of 
the small-scale field does not make a contribution to the mean emf 
£, correlated with the imposed test fields. However, small-scale dy- 
namo action makes the calculation of reliable average values of a 
and rjt more difficult. Of course, in the supercritical case the long 
time limit of any kinematic problem becomes unphysical. Never- 
theless, within a certain time interval the averaged values of a and 
rjt match reasonably with theoretical expectations of FOSA. 
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